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Abstract 



u. 

■ We consider the inverse conductivity problem with one measurement for the equa- 
tion div{{<Ti + ((T2 — cri)xD)Vu) = determining the unknown inclusion D included 
in r2. We suppose that is the unit disk of R^. With the tools of the conformal 
mappings, of elementary Fourier analysis and also the action of some quasi-conformal 
mapping on the Sobolev space H^/^(S'"'^), we show how to approximate the Dirichlet- 
to-Neumann map when the original inclusion D is a e— approximation of a disk. This 
enables us to give some uniqueness and stability results. 

> . 
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Q ■ 1 Introduction. 



In this paper, we study the inverse problem of conductivity with one measurement. Given 
a bounded domain $7 C with reasonably smooth boundary, an connected open set D 
compactly contained in fi, we consider for any / G H^/^(3Z)) the problem of recovering 
the subset D entering the Dirichlet equation 



p^j^ j^l divi{a, + ia2-a,)xD)Vu) = in J7 
- - ' L J J 1^ u = /, on oU 

from the knowledge of the current flux g = cidnU in dO, induced by the boundary value / = 
■f^l^n- We will denote Ad '■ H^/^(f?0) — > H~^/^(9r2) the Dirichlet-to-Neumann map which 
maps the Dirichlet data / onto the corresponding Neumann data g = Au{f) = aidnU. 
We can reformulate the inverse problem with one measurement as the determination of 
D from the Cauchy pair {f,g). We mention here that we do not need the full knowledge 
of the Dirichlet-to-Neumann map but only one pair of Cauchy data (/,(?)• For such a 
problem, we know that the uniqueness question is, in general, an open problem. It has 
been solved only for the special class of convex polyhedra, disks and balls. For other 
domains, Fabes, Kang and Seo have studied the global uniqueness and stability within 
the class of domains which are e— perturbations of disks. The main ingredients in the 
work were layer potential techniques and representation formula for the solution ud oi the 
problem P[D, f]. 
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Our main goal is to revisit the paper [HI with other techniques than boundary inte- 
gral representations. Throughout our paper, the two-dimensional case will be considered. 
Instead of layer potential techniques, conformal mappings and Fourier analysis will be 
another approach to review the two questions of stability and uniqueness within the class 
of disks and perturbed disks. 

Let us illustrate briefly the main steps of our arguments. Since il. is doubly connected, 
conformal mappings allows the construction of the conformal transplant function u which 
is solution of an elliptic problem that is obtained by transporting the original problem 
P[D, f] by means of a change of variables induced by the conformal mappings. A natural 
way to study the original Dirichlet-to-Neumann map is to study the transplanted Dirichlet- 
to-Neumann map; indeed one can show that when the original D is a disk in fi, then 
we can give an expression of the new Dirichlet-to-Neumann map by means of Moebius 
transforms and the classical formula of the Dirichlet-to -Neumann operator related to an 
concentric annulus. The elementary properties of Moebius transforms allow us to get 
an uniqueness result within the class of circular inclusions and for some special Dirichlet 
boundary measurement. 

When D is not a disk, things become more difficult. The conformal transplantation 
furnishes a Dirichlet-to-Neumann operator that is not very convenient to study. Indeed, 
all the classical tools of perturbation theory by compact operators and Von-Neumann 
expansion of the inverses are absent and thus we have no way to make our explicit formula 
more suitable for numerical purposes. However, when the original inclusion D is an e— 
perturbation of a disk B, then one can show that we have a reliable expression of the form 
Ad = Kb + Re with a remainder that is of order e". We show that a depends on 
the Sobolev regularity of the Dirichlet boundary measurement / and of the regularity of 
the boundary dD. In the conformal transplant, we have to deal with the two conformal 
mappings that map respectively ^\D into the annulus and and D on the ball; and the 
restriction of the maps on the corresponding boundaries will be of great importance in the 
error estimate. In our context, the diffeomorphism ^ is obtained from a composition of 
the two boundary correspondence functions. The error estimate is not straighforward, it 
is a consequence of the hardest problem of estimating — /io^||jji/2^^i^ when £^ : S"^ ^ S"^ 
is a W^'°° diffeomorphism of the circle and when /i is a function that belongs to some 
Sobolev space W{S^), s > i. We were not able to give the best Sobolev exponent s for 
which the estimate is true. However we give a result for the exponent values s = 1 + a 
for some < a < 1. At our best knowledge, the question remains open when h belongs 
to H'^(S'^) when ^ < s < 1. Our result about the precomposition of Sobolev spaces with 
quasi-regular diffeomorphisms are in the continuation of the pioneering works (see 11] ) 
where are studied the action of quasi-regular homeomorphisms on the critical Sobolev 
space Hi/2(5i). 

Let us point out that the resolution of an inverse boundary value problem for harmonic 
functions arising in electrostatic imaging through conformal mapping techniques has been 
introduced by Kress and his collaborators. The interested reader can consult the seminal 
work of Kress and al ^ Ej and our paper in 0. 

The paper is organized as follows. In section 2, after introducing some definitions 
and recalling some preliminary results concerning Moebius conformal mapping, we state 
the uniqueness results for disks. In section 3, we investigate the continuity properties 
of the superposition operators on Y{^/'^{S^) generated by regular diffeomorphisms of the 
circle. We then describe the approximation of the Dirichlet-to-Neumann map obtained 
after a sufficiently small deformation of a disk. In section 4, we prove the main result 
of uniqueness for disk and the e identifiability of e disks. In section 5, we prove the 



precomposition inequality. 



2 Main assumptions and results 

We shall assume throughout that Q is the unit ball of R^. Let us introduce the notion of 
small perturbation of disks. Given e > 0, a domain D is called an e— perturbation of 
a disk if there exists 6 G C^{dB) with ||(5||c2(g5) < 1 such that 

dD : X + e5{x)i'{x), x e dB 

where is the outward unit normal to dB at x. Denoting f^o C ^2 the set of points at 
some distance 5o from dQ, we will denote by C[e] the class of e— perturbations of all disks 
contained in Qq with the radius larger than a fixed number po that can be arbitrary small 
provided than it remains big with respect to e. 

We will assume that the domain D entering in equation is a disk or an e— perturbation 
of a disk B G i^o- Our main results concern the identif lability (the case of a perfect disk) 
and the approximate identifiability. 

In a first time, we deal with the perfect case where e = 0; we have 

Theorem 2.1 Let Di be a disk centered at the origin and of radius Ri. Then there 
exists a boundary Dirichlet measurement f{9) = cos 6 such that if D2 is an arbitrary disk 
contained in Q and A£)^(/) = h.o^{f) then Di = D2. 

In a second time, we consider the case of perturbed disks. We take the same boundary 
measurment f{9) = cos 9. We then have the following 

Theorem 2.2 Let Dq £ C[e] and let A/)^ (/) = g. Then there exists a positive constant 
C > such that if D £ C[e] and Au{f) = g on dfl, then 

\DADo\ < Ce", (2) 

where < a < 1 is a constant depending only on the a priori data and where DADq 
denotes the symmetric difference of D and Dq . 

The Sobolev 11^/^(5^) plays an important role in the proof of stability. We recall that 
H^/^(S'^) stands for the Hilbert space of real functions / defined on (modulo the con- 
stants) whose Fourier expansion 

n=oo 

fie'')= c„(/)e-^ 

n=— 00 

where the Fourier coefficients (c„(/)) are such that the sequence {\/nCn{f))n is square 
summable. For each / belonging to the Sobolev space 11^/^(5"^), its norm is the weighted 

P norm {J2'^=-oo I'^l W^e will also write that / belongs to the half Sobolev 

1/2 

space if and only if the sequence of its Fourier transform {cn{f )) belongs to I2 (Z). 

Let us recall some results about the action or composition (it is en fact a "precom- 
position") by quasi-symmetric homeomorphisms of the circle S^. Given an orientation 
preserving homeomorphism £^ : S"^ ^ oi the circle, we consider the superposition oper- 
ator generated by ^ defined by 



Fi{f) = foC, /GHV2(5i). 



It is known that maps ^^^"^{3^) onto itself if and only if is quasi-symmetric in the 
sense that we must have the doubling condition 

— - < K, 

m)\ - ' 

where X > is positive, where / is any interval on of length less than vr and where 21 
is the interval of after doubling I but by keeping the same midpoint. Furthermore, we 
have 

ll-^«ll£(Hi/2(5i),Hi/2(si)) < \l K + 

We recall also that among all quasi-symmetric homeomorphisms of , the Moebius trans- 
formations of act unitarily on }1^^'^{S^). 

An important question arises : can we hope to bound the error norm 1^1/2(51 ) 

by C||n||jji/2(5i'). The answer is important since it will allow us to estimate the error 
between the original Dirichlet-to-Neumann operator Ad and the transformed Ab- A first 
idea is to guess that such an estimate can be possible if ^ is not far from a Moebius 
transform of S^. However, at this stage of our work, we have to add some regularity 
assumptions on the target function u. To be more precise, we are only able to prove the 
following result of continuity for the precomposition by diffeomorphisms. 

Theorem 2.3 Let < 6 < 1 and let u be a function belonging to H-'^+'^. Let (f> be a regular 
quasi-regular function that we suppose to be a W^'°^ difjeomorphism on . Then we have 

||m o - m||jji/2(51) < C{6') ||u||Hi+i(5i) - ^llwi'°°(5i)) (3) 

holds for all 6' G (0, 6) where tus' is the modulus of continuity defined by 

=max(/+i/2,t'5'). (4) 



3 Change of variables and superposition operators. 

We suppose that D G C[e] is an e— perturbation of a disk. We aim to approximate 
the Dirichlet-to-Neumann operator A/j by A^. In a first time, we wish to transport the 
original problem P[D, f] by means of conformal transforms and to study the corresponding 
Dirichlet-to-Neumann operator. 



3.1 Change of variables and analysis of the Dirichlet-to-Neumann map 

Thanks to the classical mapping theorems, we know that there exists p £ (0, 1) and an 
analytic function that maps bijectively ^}\D onto the annulus 0, \ Bp where Bp is the 
centered disk of radius p. If the outer boundaries correspond to each other and if the image 
of one point on d^l is prescribed then is uniquely determined. Furthermore, thanks to 
the Riemann's mapping theorem, we know that there exists also a conformal mapping <I>* 
that maps bijectively D onto Bp. We recall that the restrictions of the conformal maps to 
the inclusion have the same regularity than dD (see jSlE] for more details). 

We denote by = ($*)~i (respectively {^'^)~^) the inverse of (respectively ^'^) 
and by 7 : [0, \dD\] — > dD the parametrization of dD in terms of arc-length. We set 



(5) 



and 

r{e)=^-U^'ipe^'')]. (6) 



Let U^^ (respectively U^^) denote the conformal transplant of iuD)\Q\~D (respectively 
(ud)\d)- We have 

-Op 

and 

= Uh^ o 

let us give the explicit form of the elliptic equations satisfied each of the conformal trans- 
plants. We have 

Proposition 3.1 Let ^ be the diffeomorphism on dBp defined by ^ = {(j)'^)^^ o (j)^. Then, 
we have 



Op 


= m \ 


Dp 


= / on (90, 


AUh 

Op 


= in Bp, 


Dp 


= U%^o^ on dBp, 




= a2drUg^ on dBp 



Proof of Proposition 13. H It is elementary and essentially based on the Cauchy- 
Riemann identities. We left the details to the reader. ■ 

3.2 Boundary correspondance for e perturbations of disks. 

The important fact to notice now is that if D E C[e\, then ^ is a perturbation of the 
identity. More precisely, one has the following result 

Proposition 3.2 There exists a positive constant C > such that 

U-I\\wi^is-)<Ce (7) 

holds for all D eC[e\. 

Proof of Proposition 13. 2t Since ^ = {<j)^)~^ o (p^, the proof is split into two parts: 
in a first time we estimate the contribution of the interior then in a second time the 
contribution of the exterior. We claim that 

W-I\\wUsn<Ce, (8) 

and 

U'-I\\w^iS^)<Ce. (9) 

Deducing (O from the claims ((SJ,® is easy and left to the reader. We now prove the 
claims in the two next sections. ■ 



3.2.1 Proof of claims dHl)-®- 

The simply connected case: claim Q. Without loss of generality, one can assume 
dio to be starlike with respect to the origin. We use polar coordinates to write 

dD: z = z(0) = r(0)e^'^, < </> < 27r, (10) 

where r is a given positive regular function of period 2tt such that 

r{e) = {R- 5{e)), < < 27r (11) 

5 being a function of period 27r satisfying 

sup <e, A; = 0,1, 2. 

0<6»<27r 

From Henrici ([H],|ni)) we learn that 9 and (p'^{6) are related by the Theodersen's integral 
equation 

cj^'{9)-9 = n{logr{m)) (12) 
where Ti., the Hilbert transform on the circle, is defined by 

1 r2iv g J. 

nf{9) = —P.V.J^ f (t) cot ^dt. (13) 

The same author learns us that the Theordersen's integral equation admits exactly one 
continuous solution (/>(^), < 9 <2tt under the condition that the ratio 



S = sup 

O<0<27r 



r'{9) 



i9) 



(14) 



satisfies 5 < 1. This condition ()14() means that the angle between the outward normal 
and the radius vectors does not exceed arctaniJ < 6. It also means that we have to deal 
with curves div that are not too far from a circle. It is referred as the 6 condition. In 
our context, D is not to far from a disk. Some straightforward arguments (primarily due 
to Montel and Lindelof) show that the boundary correspondence between 9 and (j)'^{9) is 
given by 

(j)'{9) =9-n6{9) + 0{e^). (15) 

The interested reader will find a geometric proof in ((Sj)- Hence, if the perturbation 6 is 
in W^, one can show easily that there exists a positive constant C > such that 

\(l)\9) -9\< Ce. (16) 



The doubly connected case: claim ®. As we did in the previous paragraph, we 
give the asymptotic behavior oi (j)^ — I when e — > 0. The analog of the Theodersen's 
equations for the doubly connected case is described by the so called Theodersen's and 
Garrick equations. The boundary correspondence is given by the following result. 

Theorem 3.3 Let O he a doubly connected region conformally equivalent to the annulus 
p < \w\ < 1. We suppose O bounded by piecewise analytic curves Tq and Ti, both starlike 
with respect to the origin and parametrized as above. If we suppose that 



r27r r2iv 

/ {MO)-o) d9= / iMO)-o) d9 = 

Jo Jo 



(17) 



then there holds 



-IC.ilogniMO))) 
-HpilogniMO))) 



(18) 



where the operators ICp : L'^{{0,2tt)) ^ L'^{{0,2tt)) and Hp : L'^{{0,2tt)) ^ L'^{{0,2tt)) are 
defined as follows 



and 



Amd 



Amd 



Hpie'"'') 



ICpie 



imd\ 



, m / 



m / 



Furthermore, the radius p is explicitly given by 

f27r 



(19) 



(20) 



(21) 



It is straightforward to show that the Theodersen and Garrick equations enables us to get 
@ when the radius perturbation 6(9) belongs to W'^'°°. 

3.3 Approximation of the Dirichlet-to-Neumann map for perturbed disks. 

Let D be an fixed e— perturbation of a disk B. We want to estimate the correction term 
||A-z)(/) — Ab(/)||jj-i/2(51-) with respect to the perturbation factor e. We have 

Theorem 3.4 Let D be a regular domain that is an e perturbation of a disk centered at 
the origin. Let a be a real in (0, 1). Then there exists a constant C > such that 



|A^(/)-Ab(/)||h- 



1/2(51) 



Hi+Q(5i) 



(22) 



holds for the boundary measurement f belonging to the Sobolev space lV-^°'{S^] 



The proof is lengthy and requires some preliminary results that we will state and prove 
before. 



Computation of the transplanted Dirichlet-to-Neumann map. Our main task is 
to give the analytic expression of : HV2(517) ^ R-y^{dn) defined by 



A^^(F^) = ai9.C/|^, 

where we set = f o ^Jg^ = f o tp^. The work will be divided in two parts: in the first 
one, we give some preliminary results based on the expression of h = {U%^)\QBp- This 
will allow us to get A^ {F^) and a convenient approximation Aspif)- A straightforward 
calculation shows that for p < r < 1 we have 



U'{re'^) 



Inr 



Inp 



co(M-co(n)+co(n + j:^— ^ 



|n| 



,\n\ 



in8 . 



•e\ind 



r,2\n\ 



.|n| 



c„(/i)e*"^ 



hence after identification of the Fourier coefficients, it comes that 

^ ^ ,P (23) 

^B,{F')) =ai^4J^[(l + p2H)c4F<^)-2pWc„(/i)], n^O. 

We see that the knowledge of h determines uniquely the Dirichlet-to-Neumann operator 
h-Bp- It is then useful to get some informations about h. First, we have 



Proposition 3.5 We have 



\n\cn{h o Oe*"' + 2e'(^)- -^|n|c„(/i)e 



(72 ^ 1 - p2|"l 



(24) 



Proof of Proposition 1^31 Since C/* is solution of a Dirichlet problem in the disk Bp, 
we obtain 

and thanks to the jump condition satisfied by the normal derivatives 5rf7jg^^ and drU^g^^, 
we deduce equation (|24|) . ■ 

To solve H24p is a difficult task since the explicit expression of h is hard to manipulate; how- 
ever when the perturbation factor e is very small, one can give a suitable approximation. 
Before entering in the details, we have to give some qualitative properties of h. 

A regularity result. We first prove the following tangential regularity result for solution 
of the conductivity problem. 

Lemma 3.6 Let u be the solution of the problem P[D,f] where dD is assumed to be 
of class C^. Then u^qd belongs to H^+''(9D,R) for some S G (0, 1). 

Proof of Lemma 13.61 From classical methods, the problem P[D, f] has a unique 
solution in the variational space H^(r2) with a trace u^qq = / G H^/2(30) and a normal 
derivative dnU^dn '■= 9 ^ H~^/2(9f2). For all x G dD, we use the classical representation 
formulae for harmonic functions with the help of the single layer and double layer potential: 
since u is harmonic in i}\ D and in D we have 



in+(x) = / dnG{x,y)f{y)ds{y) - I dnG{x,y)u^ {y)ds{y) 
on JdD 



2 



G{x,y)g{y)ds{y) + / G{x,y)dnU+ {y)ds{y), 
dn JdD 

u~{x) = / dnG{x,y)u~{y)ds{y) - G{x,y)dnU~ {y)ds{y) 
JdD JdD 



2 



where G is the Newtonian potential and where the normal n to dD is oriented to the 
exterior. Using the jump conditions [u] = [(T9„n] = across the interface dD, we check 
that u = = u~ solves the integral equation 



1 , , fT2 - CTl 

-u{x) + — — — 

1 (Tl+(T2 JdD 



dnG{x,y)u{y)ds{y) 



CTl + CT2 



dnG{x,y)f{y)ds{y)- / G{x,y)g{y)ds{y) 
on Jdn 



Since dD n dO, = 0, the right hand side of this equation is of class C°°. From the fact 
that the boundary dD is of class C"^, the double layer potential on dD maps }i^{dD) into 
H*+^(9Z)) for and hence is compact as operator from H*((9Z)) into itself. We conclude 
thanks to the Fredholm alternative. ■ 

The perturbation argument. Let : H^+°(9-Bp) H^^/^(9i?p) the operator defined 
by 

h ^ = \n\cn{h o oe-^ + e{0)^Yl 

and T : R^+°'{dBp) B'^/^idBp) the operator defined by 



^ 1 + 



h ^ T{h){e) = V |n|c„(/.)e*"^ + -E T 

(To 1 



,2|n 



■c„(/i)e 



in9 



We set DT^ = T^-T. We have 

Proposition 3.7 There exists a constant C > such that 

Wm^) -T{u)\\^-i/2(^QB^) < - I\\w^\\u\\iii+c.(^QB^) 

holds for all u belonging to H^^°(9i?p) with < a < 1. 

Proof of Proposition 13. 7t We decompose DT^{u) = Ti{u) + T2{u) + T3(n) where 

Tiium 



(25) 



|n|c„(u o ^ — u)e 



T2{u){e) 

Tsium 



rr 1 -I- n2|" 



<7l 
0-2 



We begin to estimate ||22(ti)||jj-i/2(gBp), we have 

ll^2(n)||H-i/2(aB^) < C(o-i,cr2)||C' - 1||L||5cIIh-i/2(9Bp), 
where 5^ = 5i o C + 52 ° C with 



= |n|c„(u)e* 



and 



While the estimation of \\g2 o S,\\ii-''-/'2{dBp) straightforward 

||c/2 0?|lH-i/2(9iJ,) < C{5o)\\u\\^i/2^QB^), 

the estimation of \\g2 o ^llH-i/2(aBp) is a httle bit harder. We first observe that, since u 
belongs to }i^{di}), s > 1, we can define Tlu' = {fisy where TC is the Hilbert transform 
on the circle. It then comes that 

(7i o ^ = Hu o ^ = (Hu)' o ^, 

and then that 

\\gi o eilH-i/2(aijp) < C\\{nuy o ^||H-i/2(aijp) < <^ll(W^i)'llH-i/2(aij^) 

< C\\nu\\-iii/2(^QB^) < C\\u\\-^i/2^gBp)- 

In the last inequality, we used the fact that the Hilbert transform corresponds to an 
unimodular multiplier and then is an isometry on H^/^. Gathering the estimates on (71 o ,f 
and gi o C, we get 

ll^2(u)||H-i/2(aBp) < C{ai,a2,So)\\^' - 1\\Io\\u\\hi/2(^oBp), 
It remains to estimate ||73(n)||jj_i/2(g^^-), we have 

1 II ^ 11^ f \ I ini(e) ine \ II 

We focus on the first part of the sum. We have 



Hence, 



(Tiu o + (1 - OC^u) o ^ - {TiuY o ^ - (?t:n)' 
((?^n) o e - Hu)' + (1 - ° 



^|n|c„H(e^"«(^)-e-^)||H-i/2(aB^) 

< ll(Wu) o ^ - Wn||Hl/2(g5^) + lie' - l||ooh'|lH-i/2(aB^) 

< II ^c„(Wn)(e-«W - e-^)||Hi/2(9B^) + ||C' - l||oo||n||Hi/2(aB,) 



Suppose an instant that u is a trigonometric polynomial of degree d. Thanks to the 
composition Theorem 12.31 for 6 £ (0, 1), we have 

||^c„(Wn)(e-«W -e-^)||H-i/2(as^) < ^ |c„(^)| ||e-«W - e-lH^/.(^5^) 

n^O |n|<d 

< C{6) ||u||Hi+a(3B^) a;5(||0 - /||wi.'-(si))> 



where u)25 is the modulus of continuity defined by Q). Then, there exists a constant 
C(<7i, C72, 5o) such that 

ll^3(^i)|lH-i/2(aBp) < Cu)s{\\i - Id\\Y^i,oo(^QBp))\\u\\ui+S(dBp)- 

The estimation of Ti(n) obeys to the same computation : we have by fohowing the same 
hnes 

ri(^)ii 

end this ends our proof. ■ 
Proof of Theorem 13.41 Let us return to the equations satisfied by h. Set 

n| 

and suppose that the perturbation parameter e is sufficiently small such that 

\\{DT{}T-^\\ < 1. 

It follows that 

h{9) = T-'b{F',p,^,ai,a2) 

r- 

T~H{F\ p, (71,(72) - T-i (/ + DT^T-^y^ DT^T-^h{F\ p, ai, aa). 
An easy calculation of shows that 

p\n\ 



T-^ {I + DT^T-^) ^(F^p,e,cTl,c72) 



He) = 2^^ .H C.(/)e-^ + H9), 

where 6u G C°° is such that for all s > ||(^ti||Hm(5i) < C\\S, — I\\wi^{S'^)\\f\\u'^/^ fo^' 
integers m € N. Plugging this expression of h in formula ((221) setting /i = it 
follows that 



Inp 



(26) 



where (5/i is a function such that ||(5/i||j|^i /2{s^) ^ Ce. 

We are ready now to finish the proof. First of all, we know from the Cauchy-Riemann 
identities that 

A*Bp(n = (5or)(V7 

and using the same arguments that we developed above, one can easily show that there 
exists a positive constant C > depending on po and 5q such that 

\\K%{F^) - Ad(/)||h-i/2(5) < Ce"ll/llHi+". (27) 



Hence we get 

||Ab(/) - Az5(/)||h-i/2(51) = I|Ab(/) - A*s,(n + A*s//) - AD{f)\\n-^/2^s^), 

< ||AB(/)-Aij/n + C7e-||/||Hi+.(si), 

< ||Ab(/) - Ab,(/)||h-i/2(si) + C|| {n{6h)y ||h-i/2(si), 

+ C'e"ll/llHi+<»(5i); 

< ||Ab(/) - Ab,(/)||h-i/2(51) + Ce"||/||Hi+.(si). 
Recall that the Fourier coefficients of Aspif) are given by 

1 + uo^l"! 

1 — ^p^\^\ 

Hence after denoting pi the radius of the disk we get 

A^(/) - Ab,(/) = ^1 ^ |n| f - ^^^l c„(/)e-^ (28) 

and this implies that we can find a constant C > depending on (5o,po ^-nd on the con- 
ductivities (Tj, 2 = 1,2 such that 

||Ab(/) - Ab,(/)||h-i/2(51) < - PiIII/IIhV2(51); (29) 
we conclude thanks to the fact that \pi — p2\ < Ce. This ends the proof of our theorem. ■ 



4 Proof of the main theorems 

We subdivide the section in two parts : in the first one, we focus on the case where the 
inclusions are disks. In the second part, the inclusions belong to C[e] . 

4.1 Identifiability for disks : proof of Theorem 12.11 

For the case where the domains are disks, we can always suppose that Di is centered at 
the origin. 

Proof of Theorem 12. H We use conformal mappings that maps a non concentric 
disk D2 into a disk centered at the origin. We know that this can be done by means of 
the Moebius transform 

z-b 

w{z) = — with \b\ < 1. 

I — bz 

In [3 El El) the interested reader will find all the details about the properties on such 
transforms . The radius of the transformed disk is denoted by i?2- If z = re*^, then w{z) 
writes p{r,9)e^^^^'^^ where 

2 r2 + |&|2-r(6e^^ + 6e-*^) 

^ '~ l + |6|2-r(6e^« + 6e-^^)' 

and where 

r sin 61 - Qbr^ + r[sme{Qb^ - m^) + 2TO6cos6l] 



(r, 6) = arctan 



2 cos 6* - 5R6r2 -m + r[(5R62 _ cj^2) ^os 9 - 25R696sin( 



Here and denote the real and imaginary part of the complex b. A straightforward 
computation shows that 



drp{l,0) 



and 5^0(1,^) = 0. 



1 + |6|2- (6e*« + 6e-i«) 
From the chain rule of differentiation, we get: 

where A^^ is the Dirichlet-to-Neumann map for the transformed and concentric problem. 
It follows that 



Adi(/) = Ad,(/)^Adi(/) 



l + |6|2-(6e^^ + 6-i«) ^ 



and the key for solving the problem is to use the Dirichlet-to-Neumann map for concentric 
disk. We have: 



l + pRl 
or equivalently: 



cos( 



1 + ^^Rf\ 



?2|fc| 



■Cfe(/' 



l + uRl 
l-pRl 



1 



(1 + \bY) cos{e) - - -{be'"" + be-'''') 



Replacing 9 hj cj) ^{0), we then get that Aoiif) = ^D-^if) implies that the 27r periodic 
function 

F{e) = (1 - E 1^1 ^^^^^(/ ° ' ^ ''''' 



l-fxRf 
1 ,v 



(1 + cos((/)-^(e)) - 3fi(6) - ^(6e2^'^ + be-^''>'~'^^^] 
satisfies Ck{F) = for all k eZ. 

We tackle the computation of these Fourier coefficients. First of all, we need to compute 
Cjk(/ o and Cjk(e2*'^ ). A straightforward computation shows that: 

( 5R(6) if A; = 0, 

Ckif ° (f>~^) = I 1 

I -(l-|6p)(-6)'=-^else; 

f bifk = 0, 
Ck{e'^~') = <j (1 - |6p)(-6)'=-i if it > 0, 
else. 



and that 



Since 



JO 



1 + 6e^' 



we deduce that 



b^ + 2b{l-\b\ 



2\ ^ , n l;,|2\2 ^ 



2ie 



1 + he^^ 



^ + (1-|&IT 



[l + be 



and then 



Cfc(e 



2i</.-l(6l)^ 



62 if = 0, 

(_5)fc-2(l _ |5|2) _ 1 _ (fe + l)|5|2j if ^ > 0, 



else. 

Hence, the Fourier coefficients {ck{F))k are given by the following formulae 

r->2 _ r>2 

co(F) = Oandci(F)=/i(l-|6|2)2 2 i 



and for k > 2: 



Ck{F) 



.2|fc| 



1 -/ii?2 



2|fc| 



2\;,fe-l 



1 - 



i + i6r)6' 



2\i,fc-l 



1 + 



(_l)fc-l(l_|6|2) 



l-fiRf 



^[k-i-\b\\k + i)] b 



-rfc-l 



{Rf^-Ri){i-m 
(i-/.flf|)(i-Mi?f: 



■fc-i 



6'=-i)(l + |6|^) 



1 



Let us show that Ck{F) = 0, V/c > implies Di = D2. First of all, the condition ci{F) = 
implies that Ri = R2- It remains to show that the Moebius transform is, in fact, the 
identity. Equivalently, we need to prove that 6 = 0. 

Let us assume that we have 6/0; the condition Ck{F) = 0, Mk >2 would imply that 
6/6 is real and then that 6 = ±6. From the identity 



2/i- 



(i?fl)(l-|6|2) 



{l-^,Rf\){l-^,Rl) 

we would get for A; = 2p + 1, p>\ 



k-l 

b\ 



(i + i&r) 



2 1 + ^lRl 



1 - 



0, 



2jA^;^^^^^ii^^^^ 

and then |6| = 1; this is impossible since the Moebius transform requires |6| < 1. Hence 
6 = is the only possibility. Gathering the two identities Ri = R2 and 6 = 0, it then 
comes that Di = D2. ■ 



4.2 Proof of theorem 12.21 : the approximate identifiabihty for approxi- 
mate disks 

In this section, we need an intermediary lemma showing the relations between the norm 
of the superposition operators generated by a diff'eomorphism ^ and its inverse Its 
proof can found in (31); we have 



Lemma 4.1 LetT^ : Ii^^'^{S^) — > H"'^/^(5"'^) be the composition operator defined by T^{f) 
f o^. Then, if is bounded on H^/^(5^), then so is T^-i- Furthermore, we have 



'«ll£(Hi/2(5i),Hi/2(s'i) 



'C-i||£(HV2(51),hi/2(51)- 



Proof of Theorem 12. 2t If ^'i is the conformal mapping which maps 0, \ Di onto an 
annulus A(l,Rf) = ^} \ B{0, Rf), it then comes from the results of the preceding section 
that 

IIAfif (/ o ^i,e) - Adi(/)||h-i/2 < Ce", < a < 1 

where (pi^s = \dn ^^'^ where Bf denotes the disk i?(0, 

Concerning D2 which is the e-perturbation of a non concentric disk B2, we begin 
to transform it via ^'2 the Moebius transform that maps B2 onto a concentric disk of 
radius i?2- It is obvious that D2 = ^'2(^2) is a shght perturbation of the concentric disk 
B2 = ^2(52). 

The chain rule derivative gives 



I62P 



I + I&2P- 626^^-626-^^ 



where (/92 = ^2 ^ |i9f^- Thanks to our preceding results, we know that if n\i?(0, i?!) is the 
conformal transform of 0,\D2 then 

I|Ab(o,R|)(/ ° fhe) - Ad2(/)IIh-i/2 < Ce", < a < 1 
and the assumption A£)-^{f) = A£).^{f) implies that 

1-1^21' 

l|AB(0,iJf)(/) - i^|j^|2_5^gie_^2g-i/Wi)(^ ° '^2)|Ih-V2(51) 
< I|Adi(/) - AB(o,Kf)(/)llH-l/2{Si) 



+ sup 

6lG[0,27r] 



I-I62 



1 + I62P - b2e' 



626 



\Af)^if oh) - AB(0,i?|)(/ o 02)||h-1/2{S1) 



where C = C((5o) is a constant depending only on 60 and where < a < 1. Following the 
same lines for the proof of Theorem 12.11 we show that A£)j(/) = K£)^{f) implies that 

I|G'IIh-i/2(51) < Ce" 
where G is the 27r— periodic function defined by 

G{e) =(1 - 1621') 1^1 1 (usLkM f ° ^ 



-l^lg^fc</.(6l) 



l + ^i{Rl) 



■e\2 



1 + COS 0-^626''' -^626- 



■2ie 



We set F = G o (f)^^; from Lemma l4. 11 we deduce that 



3^(62 



Ad,(/) = Az5,(/) 



I-^IIh-i/2(51) < Ge" 



or equivalently 



\k\ 



(30) 



We claim that this last inequality implies the two following inequalities 

< Ce^ and - < Ce^". 

Indeed, we write 



where for k > 2 



{-lf-\l-\b2\'Ylih''-'k 



Ck{F) = ak + (3k +lk, 



2\2, rfc-li 



(l-/x(^f)2)(l-^(i?|)2fc) 



l-/.(i?f)2' 

and where 7^ = —(3k- When = 1, we have 



c^{F)=^^{l-\b2\') 



2 2 (-R2e)^ - {Ri)'^ 



(l-/xi?f)(l-^i?i)' 



From (jnOJ, we get 



V- |cfc(/)P 
|ci(F)| < Ce" and V ' ,,, < Ce^". 



\k\>2 



As we got for the case of disks, the condition |ci(F)| < Ce" implies 

\{R{f-{R^f\{l-\b2\^f<Ce''^ 

where the constant C depends on 5q. Since (1 — |&2p)^ > 5i, we then get ji?! — ^fl < Ce° 
This means that the radii of the two disks are very close. 
Let us prove that 



|fc|>2 



\k\ 



Since 



Ck{F) = {-iM'^k 



we also get 



where we set 



2\2 



(l-|fc2|^) 



{Rlf^ - {R\f _ 

{I - n{R\mi - pi{Rir^)^ \[b2, 



k-l 



e\2 



Y.k\b2\'^'''^5k<e'-, 



k>2 



2\2 



(1-I&2|^) 



(i?!)^'^ - (i?f )2 



(l-/.(i?f)2)(l-/.(i?|)2/=) 



62 
62 



k-l 



1 + ^(i?f 



l + M^f) 

k{l-fiRl] 



(31) 



£\2 



A;(l - fi{Rl 



'eA2^ 



Let us check that 62 > 0. We argue by contradiction and assume the converse: if (^2 = 
then we would have 

^ 1 21 ^ (1 - ii{R\Y){i - ^\h2 ) k{i- f,{Rfr) ■ 

Taking the imaginary part of 62, we would get 62 = i(^2)~^ and then = (^f)^ which 

is impossible since we have {Rf — i?!! — C^^" with Rf,R2 > po- Hence, from we 
deduce 

|62|^ < Ce^" (32) 

where C > is a positive constant that depends on do and po- 
Let us sum up our conditions : 

• we have K-Rf)^ — {Re\) < Ce"; this means that 

ll^fl - I^PlI < Ce^" 
with a constant C > depending on 5q and pQ. 

• We have I62I ^ Ce", this means that the center 62 of B2 is near the origin 0; a 
straightforward calculus gives 

IS2ABII < Ce". 

Since I-D2AS2I < Ce, we then get 

|L>2A5|| < 

It then follows that 



and then that 

This ends the proof of the result. 



\D2ABf \ < 
\D2ADi\ < e". 



5 Proof of the precomposition theorem. 

The main tool for proving Theorem 12.31 is the following lemma. It provides the behaviour 
of the precomposition on the Fourier basis functions. 

Lemma 5.1 Assume cj) is a W^'°° diffeomorphism on S^. Then, for all 6 G (0, 1/2) 

||g.n<^(e) _ e™«||jj,/,(5,) < C{5) ni+2<5 _ /||wi.oo(5i)) (33) 

where UJ2S is the modulus of continuity defined by ^ . 
Proof of Lemma 15. H The intrinsic definition of the H^/^(S'^) gives 



Noting that 



in<l){9) g*"^ g^'^'/'C'^) _j_g*'^'^ I ^ < 2 e*"('^~"'^'^)(^) g*"(</'~-'^<^){'^) g™(</'— 

Hi/2(5i) < -^1 + -^2 



2 1 giJ^S gina 1 21 



we write ||e*"'^(^ — e*"'^||TTi/2roM < /i + ^2 where 



and 



1 I. I. \pin9 mo 1 2 

h = — |g-(</'-")(") _ Ip l^ ^ I dgda, 



h = — I — — -d9da. 

/51X51 



Concerning Ji, we follow an idea of Bourgain, Brezis and Mironescu (^^J and first integrate 
with respect to 6. Thanks to Parseval's formula, we get for all 5 G (0, 1) 

h = n f |e^"('^-")H - Ipda < C ||(/) - Idf^. 

To estimate I2, we introduce a non- negative parameter / and split I2 into + -^2 where 



jd^^ff sm^[i<p-Idm-icP-Id)ia)] 
' 2T^JJ\e-a\<i\ sin^ 

jr^^ff I sint[(.^-Jd)W-(0-/d)(a)] ^ 
' 27rJJ\e-a\>i[ sin^ 
Since for all (5' S (0, 2) and for small enough I, one has 



2 

I dOda. 

2 



sin^ [{(P - Id) (9) - - Id) (a)] 



2 



sm — 



<Cn^'ll0'-l|ll-^. 



One checks that for 6' £ (1, 2), there is a constant C{5') such that 
Concerning I2, we get easily for all 6" £ (0, 1) 

Summing up the estimates for Ii and /| and I^, we get the stated result ■ 

Let us notice that the exponent 1 + 25 can hardly by reduced since its main part Ii is 
deduced from the Parseval's equality. Another remark is that the constant C{6) blows up 
when (5 — > 0. We can now prove Theorem 12.31 

Proof of Theorem l2.3t In a first step, we assume that u is a trigonometric polynomial 
namely 

\k\<n 



We fix 6' in (0, 6) and set a = 5/5' — 1 > 0. Tlien we have 

|fc|<n 

< C(O^5'(ll0-/d|lwi.-(5i)) A:i/2+^'|cfc(7x)|. 

\k\<n 

Writing 

^ + 5' = - a5' + 1 + 5'{1 + a) = -^-a5' + 1 + 6, 

we get 

/ \ w ^ 

Thanks to the Fatou's property that is satisfied by the Sobolev spaces (^0]), we extend 
the result to the most general Sobolev space H^"'"'^(S'^). ■ 
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